Onset of Collective Oscillation in Chemical Turbulence under Global Feedback 
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Preceding the complete suppression of chemical turbulence by means of global feedback, a different 
universal type of transition, which is characterized by the emergence of small- amplitude collective 
oscillation with strong turbulent background, is shown to occur at much weaker feedback intensity. 
We illustrate this fact numerically in combination with a phenomenological argument based on the 
complex Ginzburg-Landau equation with global feedback. 
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I. INTRODUCTION 



Chemical turbulence in oscillatory reaction-diffusion 
systems can be completely suppressed by means of global 
feedback p], 0, ||. Theoretically, this fact was found 
in the complex Ginzburg-Landau equation [J i.e., 
the normal form of oscillatory reaction-diffusion systems 
near the supercritical Hopf bifurcation point [4]. Re- 
cent experiments on catalytic CO oxidation on Pt sur- 
face demonstrated the same fact, revealing also a variety 
of wave patterns caused by the effects of global delayed 
feedback [3, |5[ . A theoretical model for this reaction sys- 
tem reproduced similar behavior [H, Q . 

In the present paper, we show that yet another transi- 
tion of universal nature can occur at a certain feedback 
intensity which is much weaker than the critical intensity 
associated with the complete suppression of turbulence. 
The new type of transition is characterized by the emer- 
gence of small-amplitude collective oscillation out of the 
strongly turbulent medium without long-range phase co- 
herence. When the collective oscillation appeared, the 
system remains strongly turbulent, while the effective 
damping rate of the uniform mode (i.e., the mean field) 
shows a change of sign from positive to negative. Thus, 
the transition is interpreted as a consequence of a com- 
plete cancellation of the effective damping of the mean 
field with the effect of its growth produced by the global 
feedback. 

In Sec.HH we start with the complex Ginzburg-Landau 
equation (CGL) with global feedback. Then we derive 
phenomenologically a nonlinear Langevin equation gov- 
erning the mean field in the form of a noisy Stuart- 
Landau equation (SL). In order to clarify the nature of 
the transition of our concern, some numerical results for 
the one- and two- dimensional CGL will be compared 
in Sees. IIIII and IIV1 respectively, with analytical results 
obtained from the noisy SL. Concluding remarks will be 
given in the final section. 



II. LANGEVIN EQUATION FOR TURBULENT 
CGL AS AN EFFECTIVE EQUATION 



One-dimensional complex Ginzburg-Landau equation 
with global feedback is given by [J Q 

BtW = W+(1 + id) d 2 x W- (1 + ic 2 ) \W\ 2 W+fiW, (1) 



W(t) 



W (t, x) dx, 



(2) 



where W is a complex field, and L is the system size 
which is supposed to be sufficiently large. The intensity 
\x of the global feedback is generally a complex number. 
It is known, however, that a suitable tuning of the delay 
time in the feedback in the original system can control the 
phase of this parameter [l], 0] • For the sake of simplicity, 
therefore, we shall confine our present analysis to the 
case of real fi, which corresponds to the situation where 
the delay time in the feedback is fixed at a certain value 
but the feedback intensity is allowed to vary. A brief 
comment will be made on the case of complex ji in the 
final section. 

We first consider the system without feedback (// = 0) , 
i.e., the usual one-dimensional CGL 0,1- As is well 
known, uniform oscillations are linearly unstable and tur- 
bulence develops when the Benjamin- Feir instability con- 
dition 



1 + cic 2 < 



(3) 



is satisfied. In what follows, we will fix the parameters 
ci and C2 as c\ — 2.0 and ci = —2.0 so that the system 
may stay well within the turbulent regime. We confirmed 
that under this condition no collective oscillation exists, 
i.e., W is randomly fluctuating on a "microscopic" scale 
around the zero value without perceptible systematic mo- 
tion. The core of our argument developed below depends 
little on the choice of parameter values as far as the con- 
dition (|3J) is well satisfied. 

It is known that if the turbulence is sufficiently strong, 
which is actually the case under the above parameter con- 
dition, the system exhibits extensive chaos characterized 
by the property [9( 
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where Df is the Lyapunov dimension of the high- 
dimensional chaotic attractor describing the turbulence. 
Extensive chaos implies that the system can be imagined 
as composed of a large number of cells of equal size such 
that the fluctuations of some variables associated with 
the individual cells about their mean value are statisti- 
cally independent from cell to cell. Thus, the fluctuations 
of a macro- variable, i.e., a variable given by a simple sum 
of cell variables over the entire system, are expected to 
obey the central limit theorem. In particular, in the ab- 
sence of long-range order, the characteristic amplitude of 
W will scale like 1 / yL for large system size. If we repre- 
sent our continuous oscillatory medium with a long array 
of N oscillators with sufficiently small but fixed separa- 
tion between neighboring oscillators, which we actually 
do in numerical simulations to be described below, we 
expect asymptotically the property 

1 N ( 1 \ 

where Wj is the complex amplitude of the j-th oscillator 
in the array. The simple nature of the extensive variable 
W mentioned above also implies that its time-correlation 
function defined by 

C(t) = (W*(0)W(t)) = lim i f ' W*(s)W( s + t)d s 

(6) 

obeys a simple exponential decay law for large t, or 

CO) oc exp {--ft) . (7) 

The effective damping coefficient 7 is a quantity which 
could only be determined from the statistical mechanics 
of turbulent fluctuations which is still far from being es- 
tablished. The above decay law is confirmed from our 
numerical simulation. Figure [T] shows numerically calcu- 
lated C(t) from which the exponential law with 7 ~ 0.22 
is confirmed except for some initial transient. In this 
numerical simulation, and in all numerical simulations 
which follow, we applied an explicit Euler integration 
scheme with a constant time step At < 0.01 and a fixed 
grid size Ax = 0.5, adopting periodic boundary condi- 
tions. The system size N used ranges from 400 to 200000. 

We now introduce global feedback and study its ef- 
fects on the dynamics of the mean field. Let the complex 
amplitude be decomposed into Fourier series as 

00 

W= W k e iq i° x , (8) 

k— — 00 

where q k — 2irk/L (fc is an integer), and the Fourier 
amplitudes are defined by 

Wk = \ I We~ iqkX dx. (9) 
L Jo 
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FIG. 1: Numerically observed exponential decay of the time- 
correlation function of the mean field in semi-logarithmic 
scales. The effective damping coefficient 7, i.e., the mean 
tangent (the broken line) of the curve with the initial tran- 
sient excluded, is estimated to be 0.22. The system size is 
N = 200000. 

The uniform amplitude Wq, which is identical with the 
mean field W by definition, obeys the equation 

W Q = W - (1 + ic 2 ) J2 W kl W k2 W* ki+k2 + fiW . (10) 

fei ,k 2 

One may wish to obtain an equation for the mean field 
in a closed form, which would be a stochastic equation 
of the nonlinear Langevin type. However, deriving such 
an equation would be a formidable statistical mechanical 
problem, so that in the present paper we will content 
ourselves by simply assuming phcnomcnologically that 
the effective exponential decay of the mean field as is 
seen in Fig. [TJ is a result of renormalization of the linear 
coefficient in Eq. (|10p by the nonlinear mode-coupling 
term. If we can neglect non-Markovian effects, the result 
of such renormalization will generally take the form of a 
nonlinear Langevin equation 

W= [(/!- 7) + ioj]W + N{W) +/(<), (11) 

where a nonlinear term J\f(W) with unknown specific 
form has been included, and f(t) represents random force 
with vanishing mean. 

The analysis of Eq. (fTTj) developed in the following 
section is based on the simplifying assumption that fit) 
is white Gaussian with the only non-vanishing second 
moment given by 

(f(t)f(t'))=4r6(t-t'). (12) 

The Gaussian nature of the random force seems to hold 
due to the aforementioned extensive nature of W. The 
white-noise assumption also seems valid because we are 
particularly concerned with the situation near the tran- 
sition point where the characteristic time scale of W be- 
comes very long. We should also note that the damping 
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FIG. 2: Long-time averages of the mean field amplitude (r) 
as a function of the feedback intensity fi over a wide range of 
\i. The dotted line corresponds to uniform oscillations. The 
system size is N = 10000. 



coefficient 7 has been assumed to be unchanged when 
the global feedback is introduced, which is actually the 
property confirmed by our numerical experiments at least 
for real /i. Our final remark is that in the above argu- 
ment about Eq. (Ilip we did not explicitly refer to spatial 
dimension. Therefore, there seems to be no reason why 
Eq. (JTTJ) should not be applied to systems of two or higher 
dimensions. 

Equation (fTTj) tells that a transition occurs when the 
global feedback intensity \x becomes equal to the effective 
damping coefficient 7 of the mean filed. We denote this 
value of /j as f( c , or 



= 7- 



(13) 



Figure [2] shows numerically obtained long-time averages 
of the mean field amplitude, denoted as (r), as a func- 
tion of the feedback intensity fj, over a wide range of /i. 
Although not very clear from these data, there is an indi- 
cation of transition from vanishing to non-vanishing (r) 
at some value of /x much smaller than those giving rise 
to a hysteresis between turbulent and non-turbulent uni- 
form states. In fact, the whole analysis in the rest of this 
paper is devoted to finding out unambiguous evidence for 
the existence of a transition through a closer investiga- 
tion of small-/i region. 

In Fig. [3l two space-time plots of the complex ampli- 
tude modulus for fx = 0.1 (weaker feedback) and p, = 0.3 
(stronger feedback) are indistinguishable. In Fig.[¥J how- 
ever, two phase portraits in the complex amplitude plane 
each obtained for fi — 0.1 and /i = 0.3 are contrasted with 
each other. While the distribution of the representative 
points of the local oscillators looks almost isotropic for 
the case of /x = 0.1, implying the absence of collective os- 
cillations, such symmetry is obviously lost when \i = 0.3, 
implying the existence of collective oscillations. Qual- 
itative difference between the two situations is further 
confirmed from Fig. [5] where a trajectory of the mean 




FIG. 3: Space (horizontal) - time (vertical) plot of the com- 
plex amplitude modulus \W\ for [i — 0.1 (a) and [i — 0.3 (b). 
The system size is iV = 400. 



field over a long time at each value of fi is displayed. It 
is clear that when the feedback is weak the mean field 
is non-oscillatory, simply fluctuating (presumably due to 
the finite-size effects described above) around the origin, 
whereas for stronger feedback the same quantity clearly 
exhibits a closed orbit with some amplitude fluctuation 
again due to the finite-size effects. Thus, if there is a 
transition somewhere between these zx-values, it is pre- 
sumably characterized by a noisy Hopf bifurcation. 

The nonlinear Langevin equation (jlip actually predicts 
the occurrence of a noisy Hopf bifurcation. From the 
symmetry of our system and the smallness in amplitude 
of the collective oscillation near its onset, the nonlinear 
effects JV(W) in Eq. (fTTj) will be dominated by a cubic 
term. Then the mean field obeys a noisy Stuart-Landau 
equation 

W =[(fi- fi c ) + iu] W - (r) + ia) \W\ 2 W + f(t), (14) 

where the coefficients ix c , u>, r)(> 0), and a are real and 
depend generally on c\ and ci. Since the above equa- 
tion can be handled analytically, it would be interesting 
to compare some of the results from its analysis with 
our direct numerical simulations on Eq. (UJ) and its two- 
dimensional extension, which are the subjects of Sees. IIIIl 
and IIV1 respectively. 



III. PREDICTED CRITICAL BEHAVIOR AND 
COMPARISON WITH NUMERICAL RESULTS 

The nonlinear Langevin equation (|14[) is equivalent 
with the Fokker-Planck equation of the following form 

MM- 



J^M.*) = ;|:U-(m-*. 



— ui + ar 



d Yd 2 ' 
2 ) — P H P 

d(f> r 2 dej) 2 



nr 4 ] P + Tr-^-p\ 
or J 

(15) 



Here P(r, 4>, t) is the probability density for the amplitude 
r and the phase (f> of the mean field at time t. The above 
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FIG. 4: Phase portraits for /i = 0.1 (a) and /i = 0.3 (b) at a 
given time. The system size is TV = 40000. 



equation admits a stationary solution independent of 
given by 



Pat (r) = exp 
Various moments of r defined by 



V_ 2 M - Mc 
4r ' 



V 



(16) 
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(17) 



P st (r) rdr 



and those of the fluctuation 5r = r — (r) can be calcu- 
lated [12j |. In particular, in the limit of weak random 
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FIG. 5: Trajectories of the mean held in the complex plane 
for /i = 0.1 and /j, — 0.3. The mean field is fluctuating around 
the origin when fi = 0.1, while it forms a closed orbit with 
small amplitude fluctuations when n — 0.3. The fluctuations 
come from the finiteness of the system size. The system size 
is N = 40000. 



force, the mean field amplitude (r) near the transition is 
found to depend on /i — /i c as 




,1/2 



(M > Mc) • 

< Mc) : 



(18) 



where A is a constant. Similarly, the mean field fluctua- 
tion ((<5r) 2 ) is given by 



((Srf) 



(19) 



which holds for /i ^ /i c . It is clear that the critical expo- 
nents associated with (r) and ((5r) 2 ) obey the classical 
law or the mean field theory, reflecting the fact that the 
transition is caused by the applied mean field and not 
by the development of local order to a macroscopic scale. 
Still the transition is in some sense statistical in nature 
unlike bifurcations in deterministic dynamical systems. 
This is because the loss of long-range order is solely due 
to turbulent fluctuations, so that a full theoretical under- 
standing of the transition phenomenon would be impossi- 
ble without statistical mechanics of chemical turbulence. 

Long-time averages of r and (Sr) 2 were obtained as a 
function of the feedback intensity from numerical sim- 
ulation of Eqs. (|T|) and © with N = 200000, and the 
results are shown in Figs. [Ha) and (b), respectively. By 
assuming that the long-time averages are identical with 
ensemble averages, these numerical data were fitted with 
the theoretical curves given by Eqs. (fl8|) and (fl~9| . respec- 
tively, where a scale factor is used as the only adjustable 
parameter. Note that /j, c is not an adjustable parameter, 
but is a constant given by Eq. (jTc 
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FIG. 6: (r) 2 vs. fi (a). TV" 1 ((<5r) 2 ) -1 vs. ^ (b). The open 
circles and the solid lines are the numerical data and the the- 
oretical curves, respectively. In each of (a) and (b), a suitable 
scale factor is used as an adjustable parameter to achieve a 
good fitting between the theory and numerical experiments. 
The system size is N = 200000. 



In obtaining Eq. (fT5|) for the mean field amplitude, 
we considered the limit of weak random force. Let our 
argument be generalized to include the dependence of (r) 
on the noise intensity as well as on /j, — fi c . Because the 
origin of the random force driving the mean field is the 
finiteness of the system size, the noise intensity should 
be inversely proportional to the system size, or 



D 

N' 



(20) 



where I? is a constant independent of N. One may thus 
write the stationary distribution in the form 



P s t(r) = cxp 



4D 



(21) 



where e = [i — fj, c . Applying the finite-size scaling law 
developed in Ref. [13[ to the average amplitude of the 
mean field, we obtain a scaling form 



(r) = N 



-1/4 



F 



(eN 1 ' 2 



(22) 
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FIG. 7: (r) vs. /i for some different values of N (a). Rescaling 
of the data in (a) according to the finite-size scaling given by 
Eq. (|22p produces a universal curve independent of iV (b). 



where F is a function (called the scaling function) de- 
pending on e and N only through eTV 1 / 2 . Equation (|2"2"1) 
is a generalization of Eq. (fT8|) . Numerically calculated 
(r) for various e and N confirms this scaling law. Fig- 
ure[7Ja) shows the dependence of the long-time average of 
the mean field amplitude on feedback intensity for some 
different values of N. As is seen from Fig. [TJb), all these 
data come to lie on an identical universal curve after the 
rescalings of (r) and e by TV 1 / 4 and TV 1 / 2 , respectively. 
In this way, the finite-size scaling law (|22p is confirmed, 
providing unmistakable evidence for a phase transition. 



IV. TWO-DIMENSIONAL CASE 

The one-dimensional reaction-diffusion systems which 
we have numerically studied in Sees. [II] and IIIII are not 
very realistic. Especially, surface chemical reactions, such 
as catalytic CO oxidation on Pt, occur in two dimensions. 
As mentioned in the foregoing sections, the transition 
considered here is the mean field type, so that the na- 
ture of the transition is expected to be the same as that 
in one-dimensional systems. We carried out numerical 
simulations on the two-dimensional complex Ginzburg- 
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FIG. 8: Similar to Fig. [TJb) but in two space dimensions. 
Parameter values are the same as in the one-dimensional case, 
i.e., Ci = 2.0 and C2 = —2.0, which gives fi c ~ 0.10. This value 
is used in rescaling of the numerical data. 

Landau equation with global feedback, and obtained a 
clear indication of transition, although the correspond- 
ing value of [i c is considerably smaller than that of the 
one-dimensional case under the same parameter condi- 
tion. Figure [5] summarizes numerical results for various 
e and N, exhibited in a similar manner to Fig.[7Jb), i.e., 
in the form of (r) N 1 ^ 4 vs. sN 1 / 2 for different values of N. 
These data form almost an identical curve again, which 
we take as evidence for a transition similar to that in 
one-dimensional systems. 

V. CONCLUDING REMARKS 

Preceding the transition at which the turbulence is 
completely suppressed and uniform oscillations set in, a 
different type of transition characterized by the emer- 
gence of collective oscillations was shown to exist in the 
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